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Abstract. The entanglement entropy of a pure quantum state of a bipartite system 
is defined as the von Neumann entropy of the reduced density matrix obtained by 
tracing over one of the two parts. Critical ground states of local Hamiltonians in one 
dimension have an entanglement entropy that diverges logarithmically in the subsystem 
size, with a universal coefficient that is is related to the central charge of the associated 
conformal field theory. Here I will discuss recent extensions of these ideas to a class of 
quantum critical points with dynamic critical exponent z = 2 in two space dimensions 
and to 2D systems in a topological phase. The application of these ideas to quantum 
dimer models and fractional quantum Hall states will be discussed. 
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1. Introduction 

The entanglement of macroscopic systems measures non-local correlations of a uniquely 
quantum mechanical nature. Quantum mechanical entanglement is defined in terms of 
observing the state of a (large) subsystem of a much larger system, which will is assumed 
to be in a pure quantum state. The entanglement entropy is a quantitative measure of 
these non-local quantum correlations. 

Let be a macroscopic quantum system, i.e. a system with a thermodynamically 
large number of degrees of freedom, and A be a (macroscopically large) subsystem, 
A & Q. Let B be the complement of A, such that A and B constitute a partition of 
Q: A[J B = Q with Af]B = (l}, and F their common boundary. Let us consider a pure 
quantum state |\E') (suitably normalized (^|^) = 1) of the macroscopic system defined 
on the whole region Q. For this pure state, the reduced density matrix for region A, pA, 
and region B, ps, are the Hermitian operators 

= trB|^)(^'|, pB = trA|^)(^| (LI) 

where tiA and tr^ denote the trace over the degrees of freedom on A and B, respectively. 
The reduced density matrices pa and pb obey the obvious normalization conditions: 

tr aPa = trfipB = 1 (1.2) 

Under these circumstances, th subsystems A and B are in a mixed quantum state defined 
by their reduced density matrices. 

From their definition, it follows that the density matrices are Hermitian operators 
whose eigenvalue spectra, {A^} and {A^} (respectively), are positive and lie in the range, 
< Xt>^n < 1- Here n is a complete set of labels of the spectrum of eigenstates of each 
reduced density matrix. Thus, the eigenvalues are the probabilities of finding region A 
in a given eigenstate of its reduced density matrix pA- 

With the above definitions, the von Neumann entanglement entropy Sa for 
observing region A, is given by 

Sa = -tr^ (pa In pa) = J2\^ln\^ (L3) 

n 

and the same for subsystem B. The von Neumann entropy has the important symmetry 
property: Sa = Sb, which holds provided fl is in pure state |\E'). On the other hand, if 
the two subsystems A and B are physically separate, the wave function factorizes into 
two (pure) states with support in A and B respectively: |^) = ® In this 

case, and only in this case, the entanglement entropy vanishes, Sa = Sb = 0, as each 
subsystem is in a precisely known pure state. 

The standard notions discussed above are simple to realize in the context of small 
quantum mechanical systems. Thus, for a simple system of just two spin 1/2 degrees of 
freedom in the singlet state |^) = (l/-\/2) (| t, J,) — | i, t)), the probability of measuring 
the left spin in state | t)(| i)) is 1/2 (1/2). For this trivial case the entanglement entropy 
of the left spin (part of this singlet state) is (quite obviously) 5" = In 2. 



Entanglement Scaling at Quantum Lifshitz fixed points and topological fl,uids 3 



However, in the case of macroscopic systems as well as for quantum field theories, 
the general behavior of the entanglement entropy is far less obvious and, except in 
very simple cases, has a complex behavior. It is a non-local property of the wave 
function as a whole of a system with an infinite number of degrees of freedom. There are 
several reasons for being interested in such a non-local property of the wave function. 
In principle, in a local field theory, one expects that the behavior of local operators 
should yield a complete picture of its possible behaviors. However, the wave function 
of a quantum mechanical system (including a quantum field theory) contains non-local 
behaviors which are encoded in the entanglement. Such non-local information underlies 
the degree of difficulty of computing the wave function of an extended system when 
starting from some simple ansatz usually based on simple local properties. 

The entanglement entropy in principle depends not only on the intrinsic parameters 
of the theory, encoded in the local Lagrangian, but also on the size and shape of 
the region being observed. Let us consider the problem of the scaling of quantum 
entanglement as the region that is being observed becomes macroscopically large. Let L 
be the linear size of the entire system, and the linear size of the observed subsystem 
(A), where we will be interested in the regime L ^ La ^ a, where a is a short distance 
scale (the "UV cutoff"). It has been know for quite a long time [H [2] that free massive 
field relativistic theories in -|- 1 space-time dimensions obey a so-called area law, i.e. 
a scaling of the entanglement entropy of the form Sa = const. (La/o)"^"^, which scales 
with the area of the boundary of region A, instead of a volume scaling as with the 
thermodynamic entropy. Indeed, the main motivation of the early interest in the area 
law was the formal analogy between the area law behavior of the entanglement entropy 
and the Bekenstein-Hawking entropy of a black holejH H] Sbh = fcB^/4£p, where ks 
is the Boltzmann constant, A is the area of the event horizon of the black hole and 
ip = ^/Gh/^ is the Planck length [G being Newton's universal gravitation constant). 
However, while in the General Relativity context the prefactor of the Area Law is 
universal, in the case of the entanglement entropy is not as it depends on the choice of 
the short distance cutoff a. At any rate, the area law of the scaling of the entanglement 
entropy holds for any local field theory. f 

In this paper we will discuss the problem of universality in the scaling behavior of 
the entanglement entropy. Having noted above that, in general, for a local quantum 
field theory the leading behavior is an area law with a non- universal coefficient, we 
will be interested on whether the entanglement entropy may have universal sub-leading 
terms, i. e. universal corrections to the area law. For a generic local theory universality 
is not expected unless the field theory itself is scale invariant. In other terms, the 
quantum field theory must be a quantum critical system, i.e. a quantum field theory at 
a renormalization group fixed point describing a continuous quantum phase transition. 

The best and most completely understood quantum field theories at a critical point 

f Fermi systems at finite density (Landau Fermi liquids) are known to exhibit logarithmic violations 
the area law, lnLyi.[5l[6] This behavior is due to the existence of a Fermi surface and its quantum 
fluctuations. 
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are conformal field theories (CFT) in 1 + 1 dimensions. In this case, the area law itself 
is replaced by a logarithmic dependence on the linear size of the region, 5' = | ln(L^/a), 
whose prefactor is indeed universal, and it is given in terms of the central charge c of 
the CFT [3 [HI [9l [10] . In general space-time dimension, scale-invariant quantum field 
theories are RG fixed points. At a fixed point the RG predicts the scaling behavior 
of local operators. In practice, again with the exception of CFTs in 1 + 1 dimensions, 
such fixed point theories are generally understood in terms of various perturbative RGs, 
namely the 4 — e expansion, the 2 + e expansion, and the expansion (see, e.g. 
Ref[IIl[l2]). 

In this language there is little conceptual difference between a fixed point theory 
describing a classical thermodynamic phase transition and a fixed point of a quantum 
field theory. Indeed, in the case of relativistic quantum field theories they are related 
by an analytic continuation from Minkowski to Euclidean space-times. In essence, of 
the entire conceptual construction of the theory of quantum critical systems is based on 
this formal connection. However, in general, relativistic invariance may or may not be 
present and the theory has a dynamic critical exponent z, relating space and time, plus 
a set of scaling laws that describe the effects of thermal fluctuations (in effect, a form 
of finite size scaling) (see e.g. Ref. p[3]). 

The description of scale-invariant quantum field theories is the problem of the 
scaling behavior of local observables. Little is known in general of the behavior of non- 
local observables such as the entanglement entropy in quantum field theories. J However, 
one may ask if the scaling behavior of quantum entanglement is related in any simple 
way to the fixed point theory of local observables and, in particular, what is the theory 
of the scaling of quantum entanglement. In the case of CFTs in 1 + 1 dimensions this 
problem is by now reasonably understood, and the connection between entanglement 
entropies and the structure of the CFT is known. (For a comprehensive discussion see 
Ref.[l5]).§ 

In this paper we will discuss the problem of universality in the scaling of the 
entanglement entropy in the context of a special class of quantum critical points in 2 -(- 1 
dimensions with dynamic critical exponent z = 2: the conformal quantum critical points 
of the quantum Lifshitz universality class(es) introduced in Ref.p^. These fixed-point 
theories in 2 + 1 dimensions are interesting in that they are essentially solvable. Indeed, 
as shown by Ardonne et al[T8], these theories have the special feature the the ground 
state (or "vacuum") wave function is scale invariant. In contrast, in most scale-invariant 
theories (including CFTs in 1 -|- 1-dimensions) it is the action (and the associated path- 
integral) that is scale (and conformal) invariant but the wave function scales. The 
scale (and conformal) invariance of the quantum Lifshitz wave functions (and of their 

I For a recent review of entanglement in quantum many-body systems see Ref. |14| . 
§ Logarithmic scaling of the entanglement entropy was also found at infinite disorder fixed point of 
random quantum Heisenberg chains [TBI 1 17j . These systems that are scale invariant only in the sense of 
an ensemble average, and the physical meaning of the universal coefficient of the entanglement entropy 
is presently not understood. 
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norms) is the reason for the solvabihty of these quantum critical points. In particular, 
the weight of a field configuration of a wave function of the quantum Lifshitz class has 
the same form as the Gibbs weight of a system in two-dimensional classical statistical 
mechanics which is scale invariant as it is a classical critical point. This 2D quantum - 
2D classical connection has been exploited with great success in quantum dimer models 
at the Rokhsar-Kivelson point [HI [20] and its many generalizations [2T[ [22 [ [T8l[23l[2i [ [25] . 
The quantum Lifshitz fixed points are multicritical points separating various "valence 
bond" ordered phases (that spontaneously break spatial symmetries of the underlying 
lattice model) and are also proximate to topological phases. For this reason these fixed 
point theories typically have more than one relevant operator. In a number of cases this 
leads to the transition to become first order while in others it remains continuous. [26[ [27] . 
They are also intimately related to Z2 gauge theories and, hence, their relevance for the 
description of topological (deconfined) phases. [211 EHl 122], [29] 

The conformal structure, and their locality, of the ground state wave functions 
of quantum Lifshitz fixed points allows for a direct computation of the entanglement 
entropy and to investigate its scaling properties for this class of states[30l [311 [32] . 
Indeed, the entanglement entropy of the 2D wave functions is simply related to a 
combination of free energies of the related 2D Euclidean CFT obeying specific boundary 
conditions, whose scaling with size is well known [33]. In Refs. [30[ [32] it was found that 
the entanglement entropy at a general quantum Lifshitz fixed point scales as || 



provided the boundary F of region A is smooth and has no cusps (or other singular 
curvatures). Here a, the prefactor of the "area" law (a perimeter law in this 2D case) is 
once again non-universal, and 'Jqcp is a universal finite contribution to the entanglement 
entropy determined by the conformal structure of the wave function [32] . Eq. (11.41) yields 
the first result on the universal scaling of the entanglement entropy for space dimensions 



In contrast, the scaling of the entanglement entropy at relativistic fixed points in 
space dimensions ci > 1 is much less understood and has only begun to be considered 
quite recently. A number of important first results are given in Ref . [M] for the relativistic 
0^ field theory with a global 0{N) symmetry, whose entanglement entropy was found 
to obey the same scaling found at the quantum Lifshitz fixed points, c.f. Eq. (ll.4p . The 
structural origin of the universal coefficients in these more generic fixed point theories 
is not yet well understood. 

The other class of macroscopic physical systems for which the scaling of 
entanglement is understood are topological phases of matter in two space dimensions 
and their associated topological quantum field theories, where it is found that they obey 

II For specific geometries, i.e. a disk, there are also universal finite contributions depending on the 
aspect ratio of the chosen geometry [30l [32] . 





d>l. 
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the scaling ES] 

+ <^-^* 

where, once again, a is non-universal, and the finite term, — 7topo is universal and 
given in terms of the effective quantum dimension V = a/XIo '^a of excitations of 
the topological phase; here {da} are the quantum dimensions of excitations of type a. 
Notice the opposite sign in Eqs. fll.4p and (11.51) . 

Being topological, these theories have no dependence on any length scale {i.e. the 
correlation length is zero). Indeed, these theories are not only scale invariant but they 
have a much larger symmetry as they are effectively independent of the metric of the 
two-dimensional space on which they exist. Thus, provided the short distance structure 
is properly treated, topological phases are well described by topological quantum field 
theories. The prototype topological phases are the two-dimensional electron gases 
(2DEG) in large magnetic fields in the fractional quantum Hall (FQH) regime [37], which 
behave as topological fiuids[38l [39]. Their effective field theories are intimately related 
to Chern-Simons gauge theory [ID]. The other class of topological phases of matter (as 
yet not realized experimentally) are the deconfined phases of Z2 gauge theories and its 
generalizations [m |42l|43lll4llSlll6l|47l[24l|48]. 

We will see below that the entanglement entropy of a topological phase (and of 
a topological field theory) depends on the topology of the surface, on the state on 
the system on that surface, and on the topology of the region being observed. [49j It 
further depends on whether the state whose entanglement is being tested is one of the 
(generally degenerate) vacua of the topological phase, or if it is an excited state {i.e. 
if there are excitations carrying non-trivial labels). This feature of the entanglement 
entropy of a topological phase raises the question of whether there is a minimum number 
of entanglement entropy measurements that can uniquely determine the structure of the 
topological field theory. If this were true it would constitute a classification of topological 
field theories (at least in 2D). This is a well known open problem (see, Ref . [50l 151] ) . 

In this paper we present a review of recent results on the scaling of the entanglement 
entropy at 2D quantum Lifshitz fixed points and in topological phases and field theories. 
We begin with a brief description of topological phases and their associated quantum 
field theories (Section [2]). We then discuss the computation and scaling of quantum 
entanglement entropy for the quantum Lifshitz universality class (Section [3]). In Section 
m we compute the entanglement entropy for Chern-Simons gauge theory and for the 
associated FQH fiuids. The conclusions are presented in Section [5l 

2. Topological Phases of Matter of Topological Field Theory 

Topological phases of matter are liquid phases of electron fiuids and spin systems without 
long range order, with or without time reversal symmetry breaking. The quasiparticles 
of these fiuid phases are vortices typically with fractional charge (if the fiuid is charged 
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as in the case of the FQH states of the 2DEG) and fractional statistics[S2| which in turn 
may be Abehan and non-Abehan. 

Topological phases do not have local order parameters. Instead their observables 
are non-local objects similar to the Wilson and t'Hooft loops of gauge theories. 
This similarity is not accidental as the effective field theories of topological phases 
are topological gauge theories. Thus, they have a hidden topological order and a 
concomitant topological vacuum degeneracy which depends on the topology of the 2D 
surface on which the state is defined. [38| [53] 

The wave functions of the quasiparticles associated with topological phases are 
states that transform non-trivially under a braiding operation [M] and hence exhibit 
fractional statistics. More abstractly, the quasiparticles transform like irreducible 
representations of the Braid group. (For a good modern review see J. Preskill|55].) 
In most cases, e.g. as in the Laughlin FQH states and their generalizations, these 
representations are one-dimensional and classified in terms of a single number, the 
statistical phase. There are a number of systems, notably the non-Abelian FQH 
states [56l [571 [58| [59] in which the quasi-particles are in non-Abelian (finite-dimensional) 
representations of the Braid group. In this context, this means that the states with 
a fixed number of quasiparticles are not completely specified by giving just their 
coordinates. In fact there are a finite number of linearly independent states for a fixed 
set of quasiparticle coordinates. This defines a topologically protected finite-dimensional 
quasiparticle Hilbert space. This feature of the non-Abelian FQH states has been the 
primary motivation behind the concept of topological quantum computing^^ [611 SH] and 
it is a subject of intense research [621 [63] . 

At very low energies, i.e. at energies low compared to the quasiparticle gap, 
topological phases admit an effective field theory description: Topological Field Theory, 
e.g., Chern-Simons gauge theory, discrete gauge theory. The best known examples 
are the fractional quantum Hall fluids (described by Chern-Simons gauge theory) and 
Z2 deconfined phases, which describe quantum dimer models[22] and Kitaev's Toric 
Code^e]. 

Experimentally, the "best known" topological quantum liquids are, as we stated 
above, the FQH states of the 2DEG at large magnetic fields. Most of the FQH fluids are 
Abelian FQH states: Laughlin [37] and their generalization, the Jain states [6l]. These 
states have been know to exist as plateaus in the Hall conductance: = ue^/h where 
u is the filling fraction of the partially filled Landau level. As far back as 1998 it was 
determined experimentally, through measurements of the noise of the tunneling current 
at a constriction of the FQH fluid, that their excitations are fractionally charged. [65[ [66]. 
More recent experiments [67] have attempted to measure the fractional statistics of the 
(Abelian) quasiparticles by means of a quantum interferometer along the lines proposed 
theoretically earlier on by Chamon et al[68j. 

More recent experimental work has focused on the non-Abelian FQH states. The 
best and most promising candidate is the FQH state in the first Landau level, at filling 
fraction u = 5/2. This state has been suspected to be a Pfaffian (Moore- Read) FQH 
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state (firm candidate). There is now strong evidence for g = e/4 vortex at u = 5/2 
from shot noise at a point contact experiments [69] . and DC transport experiments also 
at a point contact [70]. There is also a well defined plateau in the Hall conductance at 
z/ = 12/5, and a possible candidate is a parafermion state (Read-Rezayi) state. f 

It has also been suggested[71j that rapidly rotating ultra-cold Bose gases may also be 
in possible non-Abelian (Pfaffian) FQH state of bosons at z/ = 1 (a relative of the Moore- 
Read fermionic state |56]) but the experiment is hard and difficult to do. The recent 
discovery [72| 173] of a superconducting state with apparently spontaneous time- reversal 
symmetry breaking in Sr2Ru04 suggesting that it is a^p^+ipy superconductor has raised 
the possibility that this superconductor may have non-Abelian half- vortices [71[ [75]. J 

2.1. Hydrodynamic Picture of Topological Phases and Gauge Theory 

We will begin with the effective field theory of the Abelian fractional quantum hall states. 
The Abelian FQH states are generalizations of the Laughlin states whose wave functions 
for a system of electrons in A^^ magnetic fluxes, at filling fraction u = N/N^ = 1/m 
(with m an odd integer) are [37] 

^^{z^,...,z^) = l[{z,-z,r e-E-l^'l'/^^' (2.1) 

i<j 

where z = x + iy are the complex coordinates of the plane and i = ^JficjeB is the 
magnetic length. 

The effective field theory of the Laughlin FQH states is an Abelian Chern-Simons 
gauge theory [76] [77] and it is usually derived through a mean field theory approach (see, 
e.g. Ref.[78]). We will follow here instead a more intuitive phenomenological approach 
of Refs. [721 [HDl [3S| based on hydrodynamic considerations. Here I follow in detail the 
approach described by Wen [39]. 

The 2DEG is a charged fluid. As such it has a conserved charge current j^, 
/i = 0, 1, 2, using a 3- vector relativistic notation. From the conservation of the current 

d^3^ = (2.2) 

it follows that in 2D the current is dual to a vector field A^, defined by 

j, = ^e.^A9^^^ (2.3) 

The vector field is a gauge field since a gauge transformation A^ A^ + d^A leaves 
the currents invariant. 

Hence, we expect to be able to construct an effective theory of the FQH fluid 
solely in terms of the vector field A^ (hereafter referred to as the hydrodynamic gauge 
field[39]). However, to construct this theory we must take into account that the FQH 
fluid is incompressible, i.e. that is a fully gapped state (except possible at the edges 

f An Abelian (Jain) 2 + 2/5 state is however a possible competitor. 

I The current experimental evidence is strong but not uncontroversial as the predicted edge states 
have not (yet) been observed. 
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of the system where its chiral edge states reside§). Since the theory is gapped and 
gauge-invariant we expect the effective action for the hydrodynamic gauge field to be 
gauge invariant (as in any gauge theory) and local. Consequently the action is a local 
gauge invariant function of A^. In addition the effective action must be odd under time 
reversal, since the 2DEG is in the presence of a large perpendicular magnetic field. 

These requirements fully determine the form of the effective action since the only 
possible form of the effective action for the hydrodynamic gauge field that is local, 
gauge invariant and odd under time-reversal is the Chern-Simons action for a gauge 
field Afj, with a U{l)m gauge invariance: 

S{A) = — d^x e,,^A^d-A^ (2.4) 

Here S is the spatial manifold and is the (compactified) time direction. Furthermore, 
the requirement that the action must also be invariant under large gauge transformation 
on a closed manifold S (such as the sphere or a torus T^) further restricts the theory 
by quantizing the parameter m to be an integerJ^O] . 

The form of the effective theory is now fully determined once the coupling to an 
external electromagnetic field is specified (consistent with invariance under ordinary 
electromagnetic gauge transformations) with a term in the local Lagrangian of the form 
£int = — ej^^^ = — ^e^ii/A^^t^''-^.'^, and that the quasiparticles (the vortices of the 
fluid) couple through their currents j"^^ to the hydrodynamic gauge field Afj, with the 
standard gauge invariant minimal coupling. This theory then predicts [39] that the Hall 
conductance is (T„, = — t-, and that its excitations are vortices with fractional charge 
q = e/m, and fractional (braid) statistics 6 = ir/m. For a system of fermions the level 
m of the Chern-Simons gauge theory must be an odd integer. 

Following Wen again[39j, we now note that a 2DEG on a manifold S with a 
boundary has chiral edge states, described by a compactified chiral boson CFT 
U{l)m with compactification radius R = l/i/m and central charge c = 1. 

With some work, the hydrodynamic description also generalizes to the non-Abelian 
FQH states [m, [59], [82]. I will not give a full description here as this will be discussed 
in upcoming publication [82j. For our present purposes it will be sufficient to note the 
structure of the effective theory in a few cases of interest. 

The wave functions for the non-Abelian Moore-Read (MR) (Pfaffian) FQH states 

are 

<If^{z,,...,z^) = Fi(-^] ll{z,~z,y e-^^\^^\'"'' (2.5) 

^ * ^ i<j 

, where Pf denotes the pfaffian of the antisymmetric matrix l/{zi — Zj). 

For the u = 1, the wave function is a MR bosonic state with q = 1. It turns 
out [m [59] that in this case the effective field theory is an SU {2)2 Chern-Simons theory. 
On the other hand, for the case of the u = b/2 fermionic pfaffian and anti-pfaffian states, 

§ For an extensive discussion see Ref. [39] . 
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the theory now has an U{1)2 charge sector and an SU{2)2/U{1) neutral sector that need 
to be consistently glued together [561 ISTlllQl [831182] . 

Furthermore, the non-Abelian states have a richer excitation spectrum than their 
Abelian counterparts, which now contains: 

• Half- vortices (denoted by cr) with charge q = e/4 (fermionic case) and q = e/2 
(bosonic case) and non-Abelian fractional (braid) statistics. 

• Vortices that are charge neutral Majorana fermions (denoted by tp). 

• Laughlin vortices with charge e/m and abelian fractional statistics vr/m 

In summary, for the non-Abelian FQH states the effective theory has an Abelian charge 
sector, a f/(l)m Chern-Simons gauge theory, and a neutral sector, which is described 
by a non-Abelian SU{2)q Chern-Simons gauge theory with a U{1) subgroup moded out 
(a "coset"). The charge and neutral sectors are glued together by the requirement that 
the states thus obtained are local with respect to the electron (as the FQH fluid is an 
electron condensate) [56]. 

Below we will use this construction to compute the entanglement entropy for 
Abelian and non-Abelian FQH states directly from the effective Chern-Simons gauge 
theory [19]. Details of the construction of these effective field theories will be given in 



2.2. Time Reversal Invariant Systems and Quantum Dimer Models 

Quantum dimer models (QDM) are time reversal invariant lattice systems, originally 
proposed to describe the quantum frustration of antiferromagnetism in a doped Mott 
insulator [19]. In that framework, the dimer degrees of freedom are valence-bond singlet 
states of nearby spins on a lattice (for a detailed description see Ref . [20] ) . Quantum 
dimer models have a special value of its couplings, known as the Rokhsar-Kivelson (RK) 
point, where the Hamiltonian can be shown have an exact ground state wave function 
has the short range "resonating valence bond" (RVB) form 



where states represented by different dimer configurations are taken to be an 
orthonormal basis. || A special, and very useful, feature of this state is that the weight 
OS a configuration, i.e. a particular covering of the lattice by dimers, is the same as the 
Gibbs weight of a classical dimer model on the same lattice, in particular, the norm of 
this wave function is equal (in this case) to the partition function of classical dimers on 
the same lattice. Using the 'quantum-classical' connection it is straightforward to show 
that the equal-time correlation functions of the quantum theory can be computed from 
the correlation functions of the corresponding classical dimer model for suitably defined 

II This condition is not satisfied for a systems with spin 1/2 degrees of freedom but it is accurate in 
SU{N) spin systems in the large N limit [28j. 



Ref. [82]. 




{C} = all dimer coverings of the lattice 



(2.6) 



{c} 



Entanglement Scaling at Quantum Lifshitz fixed points and topological fl,uids 11 



observables. Since the classical dimer model is integrable, the corresponding correlators 
are known. 

This approach has been generalized to systems with more complex degrees of 
freedom. One such generalization is a quantum version of the Baxter (or eight vertex) 
model, which has ordered and topological phases separated by lines of fixed points 
of the Lifshitz universality class[T8]. Other generalizations [25| [23] include "doped" 
interacting dimer models. There are also generalizations with non-Abelian time- 
reversal invariant states based on quantum loop and net models at their respective 
'RK' point [H mi [211 [85]. 

We will not give a detailed description of the phase diagram of QDMs here. 
However, it will suffice to say that the RK point for a system on a bipartite lattice 
describes quantum (multi) critical points (if not preempted by a first order transition), 
whose effective field theory has dynamical critical exponent z = 2 and massless 
deconfined topological excitations ('spinous' and 'visons') [T9l [22| [29 | [26l [27] . 

On the other hand, on non-bipartite lattices, the QDM at the RK point is in a 
topological Z2 deconfined phase with massive spinous and visons and a topological 4- 
fold ground state degeneracy on a torus [221 [29]. This state is closely related to Kitaev's 
Toric code [16] and to the Z2 gauge theory at its (ultra) deconfined pointy 

In the next section we present the effective field theory of the QDMs at their 
quantum critical point: the quantum Lifshitz model. 

3. Scaling of Entanglement Entropy for the Quantum Lifshitz Universality 
Class 

3.1. The Quantum Lifshitz Universality Class 

The quantum Lifshitz model is the effective field theory of QDMs at criticality |86 [ [29 ] [18] . 
A quantum (as well as classical) dimer model on a bipartite lattice admits a height 
representation in terms of a set of integer-valued degrees of freedom (the 'heights') 
residing on the sites of the dual lattice (see B. Nienhuis[87j). We will label the dual 
degrees of freedom by h. However, not all height configurations are allowed. Thus to a 
configuration on with a dimer on the x-link of the direct lattice it is assigned a height 
configuration on the dual lattice such that the heights grow by one unit in the link of 
the dual lattice does not cross the dimer on the direct lattice and drops by 4 unites if 
it does. In this fashion, the assigned configurations are in one-to-one correspondence 
to each other. Thus heights are defined mod 4. This rule will be violated if a a site 
does not have a dimer attached to it (it is 'hole') and it corresponds to a topological 
excitation in the height picture. 

It is shown in Refs. |29[ [18] that the dual height model can be coarse-grained 
resulting in a model of a scalar field if, compactified (to refiect the constraints of 
the height configurations) by the requirement that all allowed operators (including the 

% This is discussed in considerable detail in the Appendix of Ref. [31] . 
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Hamiltonian) be invariant under ip ^ (f + 27rr, where r is the compactification radius. 
The quantum hamihonian for the resulting model is 

■ ""' (3.1) 



where (p and 11 obey canonical equal-time commutation relations, [ip{x),Il{y)] = 
i6^{x — y), and /c is a (so far) arbitrary parameter. 

The term "quantum Lifshitz model" follows from the path-integral representation 
of this quantum theory in 2 -|- 1 Euclidean (imaginary time) space-time whose action is 

2 



S = d^xdr 



(3.2) 



where r is the imaginary time coordinate. This (Euclidean) action is the same as the 
free energy of a classical statistical mechanical system in 3D at a (classical) Lifshitz 
critical point between a uniform (but anisotropic) state and a modulated state |88j. 

It is trivial to show that the ground state wave function '^q[^] of the Hamiltonian 
of Eq. fl3.ip is scale invariant and given by 



d'x (V^(x))^ 

^oM oc e J (3.3) 

which has the form of a local Gibbs weight in classical statistical mechanics. Indeed, it 
corresponds to the simplest critical classical system in 2D, the free boson of Gaussian 
model. Moreover, the norm of the 2D wave function is the partition function of this 
classical critical conformally invariant system! 

f [ d'x (V^(x))2 

Furthermore, not only the norm (squared) of the wave function maps to the partition 
function of the classical critical system (wit the weight squared) but the observables 
themselves can be mapped. In particular the observables of the effective theory have 
the form of vertex operators with charges determined by the compactification condition. 
(For details and generalizations see Ref.[T8].) 

In summary, this construction is a mapping to a 2D Euclidean CFT. Under this 
mapping we find that 

• The amplitude of \(p) is the Gibbs weight of a Euclidean 2D free massless scalar 
field: scale invariant wave functions 

• At these quantum critical points the ground state wave function is scale invariant 

• The equal-time expectation value of operators in the quantum Lifshitz model are 
given by correlators of the massless free boson conformal field theory with central 
charge c = 1. 

• Time- dependent correlators: dynamical exponent z = 2. 

• By matching the correlation functions of the QDM at the RK point and Lifshitz 
models, one finds that the parameter k has to be chosen to be /c = r = 1. 
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Figure 1. The disk geometry 



• For the 2D quantum Baxter wave function [18j. one finds that k varies continuously 
as a function of the Baxter weights. 

• This construction generahzes to states with non-Abehan braid statistics [24J. 

• In general the resulting theory is a unitary Euclidean CFT. 

3.2. Entanglement Entropy and Classical Partition Functions 

Let us compute the entanglement entropy for the ground state wave function "^oif] of 
the quantum Lifshitz model, Eq. (13.31) . We will follow in detail the arguments presented 
in Ref. [30j. We will consider here the geometry of a disk (shown in Figured]) with 
Dirichlet boundary conditions at infinity. Let region A be a large disk of circumference 
£ , and B be an annular region of inner circumference i and outer circumference L, such 
that L ^ i ^ a (as before, here a is the short distance cutoff). The common boundary 
r between region A and region B is taken (for simplicity) to be a circle of circumference 



For conformal quantum critical points, the Hilbert space has an orthonormal basis 
of states |{0}) indexed by classical configurations {0}, and the ground state |\E'o) of the 
bipartite system is determined by a CFT action S: 



and expectation values in this state reproduce CFT correlators. 

To investigate the universal finite terms in the entanglement entropy at 2D 
conformal QCPs, we will rely on the approach described in the work of Fradkin and 
Moore[30]. Using the "replica trick" to compute the entanglement entropyjH [9], they 




(3.5) 



Here 




(3.6) 
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showed that the trace of the nth power of the reduced density matrix, trp^, where 
is the (normahzed) reduced density matrix of a region A, with A G B separated by the 
boundary F, for the ground state \l/o on ^4 U -B, is given by 

Here Z„ is the partition function of n copies of the equivalent 2D classical statistical 
mechanical system satisfying the constraint that their degrees of freedom are identified 
on the boundary F, and is the partition function for n decoupled systems. The 
partition functions on the r.h.s of Eq. fl3.7p are = with support on region A 

and ll^^lp with support in region B, both satisfying generalized Dirichlet {i.e. fixed) 
boundary conditions on F of A and i?, and Z^yjE = ll^olP is the norm squared for the 
full system. The entanglement entropy S is then obtained by an analytic continuation 
in n, 

S = -tr {pA \npA) = - lim |-trp:^ = - log (^^) (3.8) 

Hence, the computation of the entanglement entropy is reduced to the computation of a 
ratio of partition functions in a 2D classical statistical mechanical problem, an Euclidean 
CFT in the case of a critical wave function, each satisfying specific boundary conditions. 

In order to construct trp^, we need an expression for the matrix elements of the 
reduced density matrix {(f)^\pA\4>"^)- Since the ground state wave function is a local 
function of the field 0(x), a general matrix element of the reduced density matrix is a 
trace of the density matrix of the pure state \I/gs'[0] over the degrees of freedom of the 
"unobserved" region B, denoted by 0^(x). Hence the matrix elements of pa take the 
form 

(0-^1/3^10'^) = 

- j[D4>^]e V2 2 

(3.9) 

where the degrees of freedom satisfy the boundary condition at the common boundary 
F: 

BCr : = (f)\ = (p'^\r- (3.10) 

Proceeding with the computation of trp^, it is immediate to see that the matrix product 
requires the condition (pf = 4>'f_i for i = 1, ■ ■ ■ ,n, and 0'^ = 0^ from the trace condition. 
Hence, trp^ takes the form 



^ J BCr 



(3.11) 
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subject to the boundary condition BCr of Eg. (13.101) . Notice that the numerator, Z„ is 
the partition function on n systems whose degrees of freedom are identified in F but are 
otherwise independent. Also notice the absence of the factors of 1/2 in the exponentials 
of Eq. fimj) . 

The other important consideration is that the compactification condition requires 
that two fields that differ by 27rr be equivalent. Hence, the boundary condition of 
Eq. fl3.10l) is defined modulo 27rr. ( Equivalent ly, the proper form of the degrees of freedom 
is e**^.) This means that one can alternatively define Z„ as a partition function for n 
systems which are decoupled in the bulk but have a boundary coupling of the form (in 
the limit Xr —>■ oo, which enforces the boundary condition) 

» n 

^r = - ^ Ar5^cos(0,-0,+i). (3.12) 

Here the fields 0j extend over the entire region AU B. Thus, this problem maps onto 
a boundary CFT for a system with n "replicas" coupled only through the boundary 
condition on the closed contour F, the boundary between the A and B regions. 

For the special case of the free scalar field, one can simplify this further by taking 
linear combinations of the replica fields. Then the condition that the scalar fields (pi agree 
with each other on F can be satisfied by forming n — 1 relative coordinates (pi = (pi — 
(i = 1, . . . , n — 1) that vanish [mod 2nr) on F, and one "center of mass coordinate" field 
(p = J2^=i 'pi that is unaffected by the boundary F (reflecting the fact that nothing 
physical takes place at F). Hence, the computation of trp^ reduces to the product of 
two partition functions: 

(i) The partition function for the "center of mass" field (p; since (p does not see the 
boundary F, this is just the partition function Zaub for a single field in the entire 
system. 

(ii) The partition function for the n — 1 fields (fi which are independent from each 
other and vanish [mod 27rr on F. We denote this by [Zf^^ ^. However, the fields 
(pi on the A and B regions are effectively decoupled from each other. Hence, this 
partition function further factorizes to Z^ = Z^Z^, where Z^ and Z^ are the 
partition functions for a single field (p on A and B respectively, satisfying in each 
case Dirichlet (fixed) boundary conditions {mod27ir) at their common boundary F. 

For any CFT there exists a conformal boundary condition that generalizes the 
notion of the Dirichlet boundary condition in the free case (see below). In terms of 
the partition functions Zjy, for a field in the whole system AU B that vanishes at the 
boundary, and Zp, for a field that is free at the boundary, 

and therefore 

5=-log|^ = 
/jp 



Zp 



n-l 



log 



(3.13) 



(3.14) 
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In the last equality, the Dirichlet boundary condition at the boundary was used to split 
the partition function into contributions from A and B, each including the boundary 
with Dirichlet boundary conditions. Finally, the entanglement entropy for a general 
conformal quantum critical point is just the dimensionless free energy difference induced 
by the partition in the associated CFT[30] 

S = Fa + Fb-Faub. (3.15) 

where Fa, Fb and Faub are the free energies of the CFT associated with the ground state 
wave function, with specific boundary conditions. For the case of the quantum Lifshitz 
wave function. Equation (13.31) . Fa obeys Dirichlet boundary conditions on F, Fg obeys 
Dirichlet boundary conditions on the inner boundary F and as-yet-unspecified boundary 
conditions on its outer boundary, and Faub obeys the same boundary conditions on its 
(only outer) boundary as Fb- 

We have thus succeeded in expressing the entanglement entropy in terms of 
a combination of free energies of 2D Euclidean CFTs satisfying specified boundary 
conditions. Equation fl3.15p is actually of general validity. For the case of a 
general CFT, Dirichlet boundary conditions are replaced by fixed boundary conditions, 
which generally implies that the boundary state is in the conformal block of the 
identity [89| [90| 132] . Hence, the problem of computing the entanglement entropies of 
conformal wave functions reduces to a problem in boundary Euclidean CFT in 2D. 

To determine the scaling of the entanglement entropy with the linear size of the 
region being observed we thus need to know the same scaling for the free energies of 
Equation ( ]3.15p . The latter is a problem that has been studied for quite a long time 
and much is known about it. For a large bounded region of linear size L and smooth 
boundary, F obeys the 'Mark Kac Law'[91jt 



F = aL"^ + l3L--x\nL + 0(l), (3.16) 
6 



a general result due to Cardy and Peschel|33]. Here, a and P are non-universal constants, 
c is the central charge of the CFT, and x is the Euler characteristic of the region in 
consideration: 

X = 2- 2/1-6, (3.17) 

where h and b are the number of handles and the number of boundaries of the region. 
Hence, the logarithmically dependent term in the entanglement entropy is 

= ~- {xa + Xb - Xaub) log L (3.18) 
6 

We will see below that the 0(1) correction has a universal piece related to the "boundary 
entropy" of Affleck and Ludwig in boundary CFT [92]. 

For regions A C B the coefflcient of the log L term vanishes since in this case there 
is no net change in the Euler characteristic: 

Xa + Xb = Xaub^ = Q (3.19) 

t 'Can you hear the shape of a drum?' 
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On the other hand if under some physical process A and B become physically separate 
and have no common intersection, it follows that Xa + Xb — Xaub 7^ 0. Hence, if 
the system physically splits in two disjoint parts, then there is a logL term in the 
entanglement entropy. Logarithmic terms in the entanglement entropy also arise if 
neither A nor 5 is a subset of the other region, yet share a common boundary. In 
this case, there is a logL term whose coefficient is determined by the angles at the 
intersections. Finally, if the boundary of A is not smooth, then the coefficient depends 
on the angles ai for both regions where the boundary F is singular. (For details, see 
Ref[30]). 

3.3. Universal Finite Contributions to the Entanglement Entropy: The Quantum 
Lifshitz Case 

In this section we will consider the interesting (and generic) case in which the logarithmic 
terms are cancel out, by the argument given in the previous section. This problem was 
discussed in detail by Hsu et al[32j and we follow their treatment closely. 

If the coefficient of the logarithmic term vanishes, the 0(1) becomes universal. 
The finite term is determined by the contributions of the winding modes to the 
respective partition functions. These partition functions have been computed and 
studied extensively in the CFT literature [93l [9l]- The result depends on the topology 
of the surface and on the properties of the CFT associated with the wave function. 

For the Quantum Lifshitz universality class on a cylinder, with La,b ^ i, the 
partition function for a boson with compactification radius R on cylinder of length L 
and circumference i with Dirichlet boundary conditions on both ends, which is well 
known: [95] 



where R = \/2r'^k is the effective compactification radius (as before), and A/" is a non- 
universal regularization-dependent prefactor, responsible for the area and perimeter 
dependent terms in the free energy. (There are no logarithmic terms for a cylinder or 
a torus as their Euler characteristic x vanishes.) In Eg. (13.201) t = ij is the modular 
parameter, encoding the geometry of the cylinder, and q = e^'^*'^. The elliptic theta- 
function ^s^r) and the Dedekind eta-function r]{q) are given by 

oo oo 

^93(r)=5^g^, r/(g)=g^n(l-^")- (3-21) 

n=— oo 71=1 

The important feature of Eg. (13.201) is the factor 1/R, the contribution of the winding 
modes of the compactified boson on the cylinder with Dirichlet boundary conditions. 

Putting it all together, it is straightforward to find an expression for the 
entanglement entropy using Eg. (13. 131) . In general, the entanglement entropy depends 
on the geometry {e.g. the aspect ratios L/£) of the cylinders, encoded in ratios of theta 
and eta functions. However, in the limit La ^ i, in which the length of the cylinders 
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are long compared to their circumference, the entanglement entropy given by Eq.(??) 
and Eg. (13.201) takes a simple form 

S = i2i + \nR, (3.22) 

where /i is a non-universal constant that depending on the regularization-dependent pre- 
factor A/" of Eg. (13.201) . Hence, there is a 0{1) universal contribution to the entanglement 
entropy ■Jqcp = Ini? for the cylindrical geometry. The explicit dependence of 7qcp on 
the effective effective compactification radius R = V2kr^ shows that it is determined 
by the winding modes of the compactified boson and thus it is a universal guantity 
determined by the topology of the surface. In particular we find that the universal piece 
of the entanglement entropy, 'Jqcp, for a compactified boson is a continuous function 
of the radius R, a conseguence of the existence of an exactly marginal operator at this 
QCP. We find the similar relations for all topologies we considered. 

Hence, for the cylinder geometry the 0(1) universal term S'cyUnder eguals 

^cylinder = log -R (3.23) 

where R = ^/2kr^ is the compactification radius of the associated CFT.:j: For the RK 



Dirichlet Dirichlet Dirichlet 




'A 



Figure 2. The cylinder geometry. 

guantum dimer model, k = r = 1 and we obtain R = \/2. 

Hence, for the Quantum Lifshitz model on a cylinder, the universal finite term 
in the entanglement entropy is log ^/2. It differs from its value — In 2 in the nearby 
topological phase, which as a topological phase is in the universality class of Kitaev's 
toric code. 

The case of the case of guantum Lifshitz model on a torus can be analyzed similarly. 
We now consider the case in which the full system AU B is a torus for which the real 
part of the modulus L/i ^ 1, as shown in Figj3l The two subsystems, A and B are 
now two cylinders, of length La and Lb respectively (L = La + Lb), both with the 
same circumference i. We will thus need the partition function on a torus and on two 
cylinders (with both ends of the cylinders obeying Dirichlet boundary conditions.) The 
trace trp^ now becomes 

, n _ f zULA,i)zEALB,i) \ 

- [ z^iiLj) ) ■ (2-24) 

I Notice that the compactification radius of the associated Euclidean CFT is not the same as that of 
the 2 + 1-dimensional theory. 
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Figure 3. The torus geometry. 



The partition functions for the two cyhnders, A and B has the form of Eq. (13.201) . The 
partition function for the torus is[94l [93] 

2 



^torusV-^jt; — I ^cylinder I 2' 



{L,i)=[Z^Z^,A^,i] ) , (3.25) 



where Z^^^^^j.{^,i) is the partition function on a cyhnder of length ^ and circumference 
£, with Neumann boundary conditions at both ends: 



where t = ij and q = exp(27rir). 

In the hmit i ^ a and Lb i ^ a, the entanglement entropy for the 

toroidal geometry is 

5 = /i£ + 21nf— y (3.27) 



^ 2 , 

Hence, for the toroidal geometry, the universal term is 7qcp = 21n(fcr^) = 21n(i?^/2). 
In Eq. (13.27!) /i is, once again, a non- universal factor which depends on both the short 
distance regularization and boundary conditions (in fact, it is not equal to the constant 
we also called in the entanglement entropy for the case of the cylinder, Eq. (13.221) .) As 
was the case for the cylindrical geometry, in the case of the torus 7qcp is also determined 
by the contribution of the zero modes of the compactified boson to the partition 
functions. Thus, here too, 7qcp depends on the effective boson radius R = V2kr^. 
However, the different values of 7qcp in Eq. (13.271) and Eq. (13.221) is due to the fact that 
on the torus all three partition functions have contributions from the zero modes. 

For a disk, with the geometry of Figure [11 the finite universal term S'disk depends 
also on the aspect ratio, and the universal finite term is [32] 



-In! 

71 



+ \nR (3.28) 



We close this section with the application of these results to several problems 
of interest. As noted in previous sections several systems of interest can be mapped 
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to the quantum Lifshitz model for an appropriate choice of the parameter k and 
the compactification radius r. Systems of this type include: a) the quantum dimer 
model[T9], for which A; = r = 1, b) a generalized 'interacting' quantum dimer 
model[23l [25], for which k varies continuously but r is fixed to 1, and c) the quantum 
eight vertex model of Ref. [H] for which r = 1 and k varies continuously (long the six 
vertex line) as (here c is a Baxter weight) 



In all three cases we have mapped the problem to an Euclidean boson CFT in 2D 
with compactification radius R = y/2kr^. In all three systems, there is a perturbation 
that destabilizes the fixed line and drives the system immediately into a topological 
phase. For the dimer models, the relevant perturbation that drives the system into the 
topological phase are dimers on both sub lattices [18], whereas in the case of the quantum 
Baxter model the Baxter weights that break the continuous symmetry of the six-vertex 
model (down to a Z2) are the responsible for this transition. As shown in Ref. [31] . the 
universal term of the entanglement entropy is the same (and hence universal) in the 
entire topological phase, taking the value — In 2, associated with its value at the Kitaev 
point, the stable fixed point of the topological phase. Hence, while in the topological 
phase the universal term has a fixed (and negative) value, along the critical lines the 
universal term of the entanglement entropy varies continuously as a function of the 
parameter A;.§ Hence, in general, we should expect that the universal term of the 
entanglement entropy may jump at phase transitions. || 

3.4- Universal Finite Contributions to the Entanglement Entropy: General Case 

The results of the previous sections generalize for a general conformally invariant wave 
function defined by the Gibbs weights of a 2D (Euclidean) rational conformal field 
theory (RCFT), e.g. the Ising model, quantum loop and net models. Once again, the 
norm (squared) of the wave function becomes a partition function of an associated 2D 
Euclidean CFT. However since these are no longer free field theories the identification 
of the CFT is a little more involved [24 j. 

In this section we will compute the entanglement entropy for wave functions whose 
weights are associated with a RCFT, following the work of Ref. [32]. Here we rely on 
very basic and standard properties of CFTs which can be found in many texts, see e.g. 
Ref. [01] . A 2D RCFT has a set of primary fields have an operator product expansion 
(OPE) of the form 




(3.29) 




(3.30) 



§ This should not be a surprise since the critical lines exist due to the presence of a exactly marginal 
operator in the theory. 

II Indeed, in an ordered phase, we expect the finite universal terms to vanish. 



Entanglement Scaling at Quantum Lifshitz fixed points and topological fl,uids 21 



where the integer-valued coefficients N^^^ are known as the fusion coefficients. On the 
other hand, each primary field define a set of conformally invariant boundary conditions 
labelled by a. This defines a boundary RCFT[33]. 

A rational CFT (RCFT) is a CFT with a finite number of primaries. Under a 
modular transformation the partition function of an RCFT with boundary conditions 
specified by the action of its primary fields, transforms linearly. These modular 
transformation laws are encoded in the modular iS-matrix which is related to the fusion 
coefficients by the Verlinde formula [96] 



qi qi qb 

NL = E (3.31) 

i ^ 

The partition function for a RCFT on a cylinder of length L and circumference i, 
with boundary conditions a and b on the left and right ends respectively, Za/b, can be 
expressed in terms of the characters Xi of the RCFT: 

^./. = E^^.X,(e-^^/^), (3.32) 

j 

where iV^^ are the fusion coefficients. 

The Virasoro characters Xj ^i-re given by the trace over the descendants 1$^) of the 

highest weight state, which are obtained by acting on it with the Virasoro generators 
L_n {n > 0): 

xAe-'^"^) = e-'^/'^HrJe-^^A , (3.33) 



where c is the central charge of the CFT, Lq is the n = Virasoro generator. Here the 
modular parameter is r = i£/2L. Under a modular transformation r —^/t, which 
exchanges the Euclidean "space" and "time" dimensions of the cylinder {i.e. it fiips the 
cylinder from the "horizontal" to the "vertical" position), the characters transform as 



X^[e-^'^'')=S|x,[e-^^^/'^), (3.34) 

where Sj is the modular S-matrix of the RCFT. The modular iS-matrix and the fusion 
coefficients are related by the Verlinde formula. 

The limit of interest here is, once again, L ^ i. Under a modular transformation, 
the partition function of Eq. (13.321) becomes 

Za/b = J2^abSixj{e-^''^/^). (3.35) 

In the limit |^ — > 0, Za/b is dominated by the the descendants of the identity 1 (up to 
exponentially small corrections). Hence, in this limit, 

Za/b - E iV:, 5° xo {e-'-'n ^e^Y. ^^b (3-36) 

i i 

and In Z^/b becomes 

In Za/b = + Ins-ab, (3.37) 
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dropping UV singular (non-universal) terms. The quantity In gat in Eq. fl3.37l) is the 
boundary entropy of a boundary RCFT introduced by Affleck and Ludwig[92], where 
the "ground state degeneracy" gab is given by 

9ab = J2^abS^- (3-38) 

i 

For a cylindrical geometry, the entanglement entropy is 

S = -\og( =fii-\n( ^-^) (3.39) 

This is the main result , which shows that 7qcp is in general determined by the OPE 
coefficients N^^ (which encode the boundary conditions on the partition functions) and 
by the modular 5-matrix, S*/, of the RCFT associated with the norm squared of the 
many-body wave function at the given QCP. 

It is important to note that it is also possible to define a unitary 5-matrix that 
governs the transformation properties of the wave function itself under a modular 
transformation. This modular iS-matrix plays a central role in 2D topological phases 
and in topological field theories. [40l [35l [97] However, only for topological theories these 
two S*- matrices are the same and in general are different or even not be defined at all!) 
We will come back to this issue in the discussion section. 

A particularly simple result is obtained for the case of a cylinder with fixed boundary 
conditions on both ends. In this case, Za, Zb and Zaub are cylinders with fixed 
boundary conditions, and hence the boundary states for all three cases are in the 
conformal block of the identity 1. Since in this case the only non- vanishing OPE 
coefficient is A^qq = 1, the universal term of the entanglement entropy, 'Jqcp, depends 
only on the element 5*0 of the modular iS-matrix of the RCFT: 

-fQCP = - In Si (3.40) 

Hsu et al[32] give a detail application of these results to the computation of the 
entanglement entropy in a few interesting non-trivial systems: a) a 2D quantum system 
whose wave function has the weights of the Gibbs weights of the 2D Ising model, b) 
quantum loop models, and 3) quantum net models (including the chromatic polynomial 
model ). For details see Ref . [32] . 

4. Quantum Entanglement Entropy and Chern-Simons Gauge Theory 

We now turn to the problem of computing the entanglement entropy in topological 
phases using the effective topological field theory. Here we will focus on the case of 
Chern-Simons gauge theories and the related FQH fluids. We will follow closely the 
results of Dong et al[49J. 

In previous sections we showed that the FQH wave functions represent topological 
fluids with a finite correlation length ^ oc £ (£ is the magnetic length). Recently, the 
entanglement entropy of FQH states has be computed numerically by K. Schoutens 
et al[98l [991 1100] and by Li and Haldane |101] . Here we will show {49 J that one can 
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compute the entanglement entropy directly from the effective field theory of all FQH 
states: Chern-Simons gauge theory. This result can be applied directly to all known 
FQH states. Since it uses the effective topological field theory, it computes only 
the topological invariant piece of the entanglement entropy. The resulting universal 
topological entanglement entropy is given in terms of the modular i5-matrix of the 
effective Chern-Simons theory and of its conformal blocks. 

In Ref. |l9] we computed the entanglement entropy for a level k Chern-Simons 
theory on a smooth manifold with any number of handles, using the seminal results of 
Witten[l0] for the Chern-Simons theory. The action of a (non-Abelian) Chern-Simons 
gauge theory is 

S{A) = ^ Jti (^AAdA+'^AAAAA^ (4.1) 

where, as usual, A^ is a vector field taking values in the algebra of a (compact) gauge 
group G. Here we will be primarily interested in the vase of G = SU{2). 

We will need a few important results on the structure of the Chern-Simons theory 
and its solution. 



• Following Witten[40]. we realize the states on a closed 2D surface as a path integral 
over a 3D volume. 

• Witten showed that the Chern-Simons states on a spatial manifold S (which we 
will take to be closed) are in one-to-one correspondence with the conformal blocks 
of a Wess-Zumino- Witten (WZW) CFT. 

• He also showed that the ground state degeneracy depends on the level k and on the 
topology of the surface S. 

• The partition functions {i.e. the value of the path integral) depend on the matrix 
elements of the modular iS-matrix, e.g. the partition function on with a Wilson 
loop in representation pj is 

Z{S\p,)=Sl (4.2) 

• Here the modular iS-matrix defines how the degenerate ground states on a torus 
transform under a modular transformation 

We will need some properties of the modular S matrix and of the conformal blocks. 
For the gauge group U{l)m, n = 0, . . . , m — 1, the modular S matrix is 

5. ,,M = ^e2W/m ^4_3^ 



'm 

whereas for the gauge group SU{2)k, = 0, 1/2, . . . , k/2, the modular S matrix is 



We also need the definition of the quantum dimensions dj 
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which measure the rate of growth of the degenerate Hilbert spaces of particles labeled 
by the representation pj. 

The Chern-Simons path integral (its partition function) on various manifolds can 
be reduced to its computation on a sphere using the method of (Chern-Simons) 
surgeries [ID]. Using surgeries it is shown that if a 3-manifold M is the connected sum 
of two 3- manifolds Mi and M2 joined along an S"^, then the Chern-Simons partition 
functions on these manifolds are related by 

Z{M)Z{S^) = Z{Mi)Z{M2) (4.6) 

In particular, if M is Mi and M2 joined along n S^'s, the resulting partition function is 

, , Z(Mi)Z(M2) , , 

^(^) = ^ ^(^3). (4-7) 

We will compute the entanglement entropy for Chern-Simons theory using the path 
integral approach of Calabrese and Cardy[9], suitably adapted for the system at hand 
by Dong et al[l9]. In other terms, one uses a path- integral to compute the nth power of 
the reduced density matrix. This leads to a "foliated" 3D manifold shown in Figure H] 
(for 72 = 3) 





Figure 4. Left: the cut-manifold defining the reduced density matrix pA- Right: the 
fohated manifold used to compute Trp^. Here ^ 00 is the inverse temperature. 



Let us consider first the simplest case, in which the spatial manifold is a sphere, 
S = 5^, and hence the space-time manifold is just a 3-sphere, S x S*^ = 5*^. The Hilbert 
space on 5^ is one-dimensional. Using the method of surgeries Dong et al[39] considered 
the case of S"^ with one A — B boundary {i.e. the observed region is a hemisphere). The 
two regions A and B are two hemispheres (disks). The 3-geometry is a ball. 

To construct trp^ we glue 2n such pieces together. When glued together 5^ rotated 
about the axis which has the topology of S^. For n > 2, the S*^ is obtained by sequentially 
gluing 2n disks, we find that the (normalized) trace of is 



(Z(53))'"" 



cl— n 
'^00 



In the replica limit, n — 1, we obtain for the entanglement entropy 

In Sqo = —\nV 



(4. 



(4.9) 
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which is the well known result of Kitaev and Preskill[3S], and Levin and Wen[3B] for 
the universal topological entanglement entropy. This result also holds for surfaces 
with arbitrary topology provided the region being observed A is topologically trivial, 
regardless of the pure state labeled by the representations pj. For the case of a sphere S'^ 
and a disconnected connected region A with M boundaries we find S\ ' = M In iSqo = 
-MlnP. 

Let us compute the entanglement entropy for a Chern-Simons theory on the torus 
with more than one A — B boundary. For a torus split into two regions with more 
than one (say two) boundary, we have two cases, shown in Figure [51 




Figure 5. Two cases for the torus. 

For the trivial state (no Wilson loop) the entropy is the same in both cases, 

5^(T2,2) = 21n5oo (4.10) 

If there is a Wilson loop with a non-trivial representation pj, we obtain the same result 
for the case depicted in Figure [5] (left). But, for the case of Figure [5] (right), for a Wilson 
loop in representation p we obtain instead, 

5^(r2,2,p) = 21n5op (4.11) 

For a state which is a linear superposition, lip) = J^pi^plp)^ fi^d 

Sa{T^ 2, ^) = 2 ln5oo -J^dl 1^ ^) (4.12) 

Clearly, the entanglement entropy now depends not only on the effective quantum 
dimension T) = Sqq but also on the quantum dimension of the excitation labeled by the 
representation p, as well as on which particular linear combination of the degenerate 
vacua on the torus in which the system is prepared. 

Let us consider now the computation of the entanglement entropy in the presence of 
quasiparticles, i. e. for a manifold with punctures each carrying a specific representation 
label. Let us consider the case of four quasiparticles on S"^: with four punctures. 
We will consider first the case of the Chern-Simons gauge theory SU{N)k, with N >2 
and k > 2, with two punctures carrying fundamental a and 2 anti-fundamental a* 
representations. 

If there is only one puncture in A, we find|l9] that the entanglement entropy is 
given by 

^^ = ^5-0^ (4.13) 
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On the other hand, if there are two punctures in A we have two possibihties: 

• Case I: There is a pair of a and a* in A and in B. Each pair can fuse into the 
identity or into the adjoint. For k > 2, the Hilbert space on S*^ with 2 pairs of a 
and d*'s is two dimensional. The entanglement entropy depends on the quantum 
dimensions of the conformal block. 

• Case II: There are two d's in A and two a*'s in B. The entropy now depends on 
which channels (representation) the quasiparticles fuse and on the choice of state 
(conformal block). 

Let us consider first Case I. We begin with a pure state 



where the basis states and |02) are shown in Figure [61 For a mixed state with a 



10) 



a|0i) + 6|02) 



(4.14) 




Figure 6. Two states represented by Wilson lines connecting punctures. 



reduced density matrix 



= aa*pii + ab*pi2 + a*bp2i + 



(4.15) 



the entanglement entropy is 



Sa = In So 







AilnAi - (4- l)A2lnA2 



(4.16) 



where 



\ada + ftp 




(4.17) 




ad^ + b\^ + {dl-lW 




(4.18) 



We can use instead an orthogonal basis 
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\(f)[) and 102) are conformal blocks associated with the trivial and adjoint representation 
9, respectively, which appear in d x d*. We have just calculated the fusion 
matrix. Here, ](!)[) and l^g) are the conformal blocks in one channel, — y=\(f)2) and 

^ ' ^1^2) + c^d|0i)) are the blocks in the other channel. The fusion matrix F 



day/So°y/dl-l 

is given by 

1 ^ ^ v/(^ 



-T- \ I — — I (4-20) 

da is the quantum dimension of the fundamental representation d, and a/c?! — 1 = d^, 
is the quantum dimension of the adjoint representation 6. 

Let us now turn to Case II. We begin again with a pure state 

|0)=a|0i)+6|02) (4.21) 
with |0i) and |02) depicted in Figure [3 




Figure 7. 

The computation is similar to Case I except that now the Wilson loops are braided 
and we now fuse two fundamentals, a. The representations cut by the boundary are the 
symmetric a = 6 and antisymmetric u = rank two of two fundamentals. The entropy 
still has the same form with the quantum dimensions given by 

...IM. (4.22) 

where I use the quantum-group notation [x] = ^1/21^-1/2 ; with q = e-2'r«/(^+fc), Por 
SU{2)k the fusion matrix is the same as before. 

4.I. Entanglement Entropy of FQH states from Chern- Simons Gauge Theory 

We end the discussion of the Chern-Simons theory with an application to the 
computation of the entanglement entropy in FQH states. 
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4-1.1. Laughlin States: We begin with the case of the (Abehan) Laughhn states at 
fining fraction u = 1/m. Their effective field theory is, as we discussed earher in this 
paper, a U{l)m Abehan Chern-Simons gauge theory. The resuhs we just derived teU 
us that for Laughhn states the the quantum dimensions of aU its m quasiparticles are 
di = 1 (as they are aU Abehan). Hence, Sqq = l/^/rn, and the universal topological 
entanglement entropy is 

5'Laughiin = - In y/m (4.23) 

It is straightforward to show the the entanglement entropy for any Abelian FQH state 
is 

5'Abelian = -^Infif (4.24) 

where g is the ground state degeneracy of the Abehan FQH state on the torus. 

4.1.2. Coset [SU{2)/U{l)]k theories: Chern-Simons theory whose CFT is the coset 
[SU{2)/U{1)]2 describes the two-dimensional time reversal breaking superconductors 
with symmetry + ipy. In some sense this is the simplest system with non- Abelian 
statistics. Here we will consider the general case of the coset [SU{2)/U{l)]k. 
In Ref. Hni it is shown that the modular S matrix for this coset is 



k{k + 2) 



sm 



Tiji +!){£' + I) 
k + 2 



g-i7rrr'/fc_ ^^25) 



For the case of most physical interest, we have k = 2, and the coset primaries may 
be taken to be (0; 0), (1; 1) and (0; 2). This is in fact just the chiral Ising model with 
(0; 0) ~ /, (1; 1) ~ (T and (0; 2) ~ The 5-matrix for (suj2)/U{l)) is 



2 



O cosct 

<-'k=2 



^ ( I V2 1 \ 
V2 -V2 



(4.26) 



\ 1 -V2 1 J 

which agrees with the results of Ref. [102] . 

Using the results we presented for Chern-Simons theory, we find that if the region 
of observation is simply connected, then the entanglement entropy is = In Soo = — In 2 
(in this case). Our expressions earlier in this section then yield the entanglement entropy 
in non-trivial topologies and in the presence of quasiparticle. 



4.1.3. Moore-Read and Read-Rezayi FQH states: pfaffian and generalized parafermion 
states: We now turn to the Moore-Read and Read-Rezayi non-Abelian FQH states, 
and their generalization. The filling factor of these states is z/ = k/{Mk + 2); M 
even corresponds to bosonic states and M odd to fermionic states [58j- These states 
are described by [SU{2)/U{l)]k x U{1) CFTs, with a suitably defined level for the 
U{1). Examples of these states are the well known Moore- Read pfaffian states. The the 
fermionic state with k = 2 and M = 1 has filling factor 1/2 (5/2 in the experiment), 



Entanglement Scaling at Quantum Lifshitz fixed points and topological fluids 29 



and the related bosonic state at filling factor v = 1 has k = 2 and M = 0. The states 
with k > 2 are the Read-Rezayi parafermionic states. 

We will discuss both the general fermionic and bosonic states with fixed k and M. 
The RCFT of interest is in all cases embedded in {SU{2)/U{l))j^ x t^(l)fc(Mfc+2)- ^^^^ 
consider the cases of k even and odd as their structure is somewhat different. Here we 
only present details for the simpler cases. The details of the derivations for the general 
case are given in Ref . |49j . 

By reasoning similar to the above, the resulting iS-matrix can be obtained by 
multiplying coset and U{1) characters. For the pfaffian state k = 2, the coset is a 
Z2 parafermion. The resulting iS-matrix will, up to identifications, be given by 

'^'^ = {Sr^)^e,r/''''\S^^'^^-'^^)/. (4.27) 

Primaries of this theory will be given by products of the Z2 primaries {/,cr, ■?/'} with 
f^(l)4Af+4 primaries of the form Oi/p. We seek a set of such operators that close under 
operator products and are local with respect to a suitable extended current algebra, 
which will be generated hj J± ^ ip e^*^*^"^^"^, where ip is the Majorana fermion of Z2. 
For simplicity, we will consider two cases here, M = (take p = 2, p' = 1, radius R = 1) 
and M = 1 (take p = 4, p' = 1, radius R = \f\j2). 

In the case of = 2 and M = 0, we find the integer-weight J± ^ e^*"^ as 
suitable extended currents. Requiring locality of operator products, we then find that 
the primaries of this theory are given by /, '0, cxe*'^''^ (all others are related to these by 
action of J±). These in fact are just the primaries of SU{2)2, as we should expect. This 
is the bosonic pfaffian state. 

In the case of k = 2 and M = 1 (the fermionic pfaffian state), we find J± ^ ip e^^^'f' 
as suitable extended currents [56]. Requiring locality of operator products, we then find 
that the primaries of this theory are given by 

/,^,ae±^<^/2V2^e±^<^/^. (4.28) 

This set closes under fusion (up to the action of J±). These operators have weights^ 
(0,0; 0), (0,2;0), (1, 1; ±1) and (0,0; ±2) respectively. We can then read off the S- 
matrix: 





( 1 


1 


V2 


V2 


1 


1 \ 




1 


1 


-V2 


-V2 


1 


1 


1 


V2 


-V2 








+iV2 


-i^/2 


2V2 


V2 


-V2 








-iV2 






1 


1 


iy/2 


-iy/2 


-1 


-1 




V 1 


1 


-i^/2 


+zV2 


-1 


-1 / 



(4.29) 



from which one can read-off the total quantum dimension is P = 2\/2. 

We will now consider the interesting example of the parafermionic states at = 3 
and M = 1: the Read-Rezayi parafermionic state for fermions at filling factor 2 + 2/5. 

^ The notation (i',r; s) represent the coset weights (^,7") and the f7(l)-charge s. 
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The k = 3 coset has primaries at {i,r) = (0,0), (1,±1), (2,0), (3, ±1), which we will 
refer to as I,o'±,e,ip± respectively. Explicitly, denoting Sp = sin(7rp/5), we have 



Ocosct 

'-'fc=3 



2 

Is 





Si 


S2 


S2 


S2 




Sl 


Sl 


\ 




S2 




g+i7r/3_g_^ 


-Si 


— e" 










S2 






-Si 


— 










S2 


-Si 


-Si 


-Si 




S2 


S2 






Si 




_g+i7r/3^2 


S2 


— e" 








V 


Si 


_g+i7r/3^^ 


_g-i7r/3_g^ 


S2 


— e" 


Htt/3 „ 
^1 


_g-i7r/3_g^ 


/ 



(4.30) 



For this case there is an extended algebra generated by the h = 3/2 operator 



ipj^ g5i0/vT5^ where is a free boson of the U{1) theory that we are attaching. Rep- 
resentative primaries are (£, r;s) = (0,0; 0), (3,— 1;1), (3,1; 2), (0,0; 3), (3,-1; 4) and 
(£, r; s) = (2, 0; 0), (1, -1; 1), (1, 1; 2), (2, 0; 3), (1, -1; 4). One can check that these have 
local OPE's with and are closed under fusion. As we will see, it is convenient to 
group them into groups of + 2 = 5, as given. The theory obtained this way is actually 
an = 2 superconformal theory, with supercharges Q± (Q_ being tp_ e-5*'/'/v^^_ 
groups collections of conformal primaries together, i.e., {(0, 0; 0), (3, 1; 5), (3, —1; 10)}, 
{(3, -1; 1), (0, 0; 6), (3, 1; 11)}, {(3, 1; 2), (3, -1; 7), (0, 0; 12)}, 

{(0, 0; 3), (3, 1; 8), (3, -1; 13)}, and {(3, -1; 4), (0, 0; 9), (3, 1; 14)} and {(2, 0; 0), (1, 1; 5), (1, -1; 10)}, 

{(1, -1; 1), (2, 0; 6), (1, 1; 11)}, {(1, 1; 2), (1, -1; 7), (2, 0; 12)}, {(2, 0; 3), (1, 1; 8), (1, -1; 13)}, 

and {(1, —1; 4), (0, 0; 9), (1, 1; 14)}. Each of these triplets represents a superconformal 

family. When we compute the iS-matrix with respect to the extended symmetry, we 

treat these groupings as one. That is, computing the iS-matrix element on the grouping 

gives a 3 X 3 identity matrix times a factor. We collect those factors into the following 

iS-matrix. 



qFRR 
'->fc=3 



sin(7r/5) sin(27r/5) 
sin(27r/5) — sin(7r/5) 



/ 1 

1 
1 
1 

V 1 



1 

UJ2 
UJ4 
UJi 
UJ3 



1 1 

UJ4 UJi 

UJ3 ^2 

UJi UJ4 



1 \ 

Us 
Ui 

UJ2 ) 



(4.31) 



where we have used the ViV) 5-matrix is S^^ = _i^e27r«ss7i5_ Above we used the 
notation is tUp = e^'^*^/^. The coefficient out front is • ■ the factors being the 
coefficients of the coset 5-matrix, the U{1) 5-matrix and the order of the automorphism 
(5 in 15), respectively. Note that it is easy to read off then the total quantum dimension 

^ = ^ = 2imW = ^ = ^5 + = (4.32) 

where here (p = {^/E + l)/2 denotes the Golden Ratio. 

In Ref . |49] it is shown that for general k and M, the primaries are the highest 
weight states of the form 

+ nk-{Mk + 2)[§] 



i^ie/-2i^]) exp z- 



^k{Mk + 2) 



(4.33) 
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or 



.e + nk-{Mk + 2)[^ 



, ^ , , (4.34) 

y/k{Mk + 2) ' 

where Zfe-parafermion primaries, n and i are integers , [x] is the closest integer 

to X. For general k the <S-matrix is given by 



S{e;n} 



{e';n'} 



y/{k + 2){Mk + 2) 



sm 



nji + !){£' + ly 
k + 2 



exp 



m{-MU' + 2ln' + 2f n + 2knn') 



Mk + 2 



(4.35) 



One can read off from this the total quantum dimension V for all M and k, since 



o 

lue 

region can be obtained. 



sm 



TT 



(4.36) 



T> " V('^ + 2)(M/t + 2) + 
from where the value of the topological entanglement entropy for a simply connected 



5. Outlook 



We discussed the behavior of the entanglement entropy near quantum phase transitions 
and in topological phases. In some special cases, the entanglement entropy of 2D QCPs 
with conformally invariant wave functions has a universal logarithmic terms. However, 
if the logarithmic term is absent the 0(1) term is universal. We also showed that in a 
topological phase the finite term in the entanglement entropy is a universal property of 
the phase, which in general is different from its value at the QCP. 

We computed the topological entanglement entropy for Chern-Simons gauge 
theories. The entanglement entropy of abelian and non-abelian FQH states is given 
in terms of the modular S matrix for the effective CS theory. This requires to glue U{1) 
charge sector and the coset (S'f/(2)/f/(l))^ neutral sector. 

For a simply connected region it is universal and depends only on the total quantum 
dimension. For regions which are not simply connected, the entropy is additive. The 
entropy of disjoint regions on a torus depends on the effective quantum dimension and 
on the state on the torus. The entropy for a simply connected region on the sphere with 
4 quasiparticles (punctures) depends on the conformal block. The entropy depends on 
the quantum dimensions and on the particular state that is chosen. The change of basis 
depends on the fusion matrix and on the conformal weights as well. 

These results suggest that it may be possible to determine the structure of the 
topological field theory by means of entanglement entropy measurements. 

We close with some comments on the question on how (or even whether) the 
entanglement entropy may be measured experimentally. Our results indicate that in 
the case of the conformal quantum critical points in 2D there is an intimate relationship 
between the entanglement entropy and the Affleck-Ludwig vacuum degeneracy of the 
associated Euclidean CFT. On the other hand the work by Fendley, Fisher and Nayak 
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on entanglement at point contacts of the non-Abelian FQH states pUB] also relates the 
entanglement entropy (change) of the bulk topological fluid, to the (change) of the 
AfHeck-Ludwig entropy of the point contact, a property of the edge states. These 
seemingly unrelated results suggest that there may be a deeper connection. More 
important perhaps changes in the bulk entanglement entropy can be measured by 
monitoring properties of a suitably defined point contact. Recent work by Klich and 
Levitov [104j in the context of a simple (and hence solvable) system of free fermions, 
suggest that this may be possible. We have preliminary results that support this idea 
in a more general context [105] . 
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